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Abstract 

We analyze the Mayer and virial series (pressure as a function of 
the activity resp. the density) for a classical system of particles in 
continuous configuration space at low temperature. Particles interact 
via a finite range potential with an attractive tail. We propose physical 
interpretations of the Mayer and virial series' radii of convergence, 
valid independently of the question of phase transition: the Mayer 
radius corresponds to a fast increase from very small to finite density, 
and the virial radius corresponds to a cross-over from monatomic to 
polyatomic gas. Our results are consistent with the Lee- Yang theorem 
for lattice gases and with the continuum Widom-Rowlinson model. 

1 Introduction 

The present work started from a seeming contradiction between results on 
cluster size distributions at low temperature and low density [TU] and pre- 
dictions from the Mayer activity expansions. It turned out that not only is 
there no contradiction, but moreover the interplay between the two differ- 
ent approaches considerably helps the physical interpretation of the classical 
expansions. 

The seeming contradiction is the following. Consider a classical system 
of particles, interacting via a stable potential with an attractive tail. As is 
well-known from the theory of Mayer expansions (see e.g., the classical text- 
book [18j), at low density, the system behaves approximately like an ideal 
gas, suggesting that particles move more or less independently and are typi- 
cally far away from each other. But [10] showed that when both the density 
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and temperature are small, there can be regimes where particles form small 
compounds - the system could behave, for example, like a diatomic gas. Such 
a behavior is, in fact, well established for quantum Coulomb systems [31E1E]- 

In order to fit both pictures together, we investigate the temperature de- 
pendence of the Mayer and virial series. The temperature dependence con- 
firms the intuitive relation between Mayer's series and Frenkel-Band theory of 
physical clusters, as exposed in P Chapter 5] : the Mayer series /5p = ^ bkz'' 
looks like the pressure of an ideal mixture of size k components with re- 
spective activities bkz'^. But unlike physical activities, the coefficients bk can 
be negative. The gas is therefore, at best, an exact ideal mixture of ficti- 
tious objects, ^^mathematical" clusters. On the other hand, at low density, 
it is tempting to consider the system as an approximately ideal mixture of 
''physical" clusters, groups of particles close in configuration space. 

Each physical cluster comes with a partition function over internal degrees 
of freedom. At low temperatures, it is natural to approximate the internal 
partition function as exp(— /JE'^), with Ek a ground state energy, and we 
expect 

Pp^J^z^expi-PEk). (1) 

k 

We prove that at low temperatures, the Mayer series coefficients bk{f3) indeed 



behave as exp{—j3{Ek + o(l))) (Theorem 3.1), so that the approximation de- 
scribed above matches the exact Mayer series. As a consequence, we can 
easily understand the formation of compounds: if /3 — t- oo at fixed chemical 
potential /i, we have to maximize {kfi — E^) in order to see which k gives the 
dominant contribution. In particular, even when the Mayer series converges, 
at low temperatures the main contribution does not necessarily come from 
k = l. 

With the approximate formula ([T]) in mind, we prove several results on 
Mayer and virial series and low temperature statistical physics, stated in 
Sect. [3j the proofs are given in Sects 4-7. The main results for the pressure- 
activity expansion are the following: Theorem |3.1 relates the Mayer coeffi- 



cients bk{(3) to the ground state energies E^, justifying Eq. ([T| as described 



above. Theorem 3.5 shows that the radius of convergence of the pressure- 
activity expansion is approximately exp(/3eoo) where e^o = lim.{Ek/k), as 
expected from Eq. ([T|. Theorem 3.2 shows that the value /i = Cc 



log R^^^{/3) corresponds to a cross-over from an exponentially small den- 
sity to a positive density, thus giving a soft physical interpretation to the 
radius of convergence. 
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In the same spirit, Theorems 3^ and |3.8 give an interpretation of the ra- 
dius of convergence of the pressure-density series. Here there are two possible 
scenarios: either the monatomic gas condenses right away to a sohd, or there 
is an intermediate phase of a polyatomic gas. In the latter case the radius 
of convergence of the virial series is given by the cross-over from monatomic 
to polyatomic gas; in particular, the virial series ceases to converge before 
the Mayer series does. Props. 3.10 and 3.11 describe the low-temperature 
asymptotics of the coefficients of the virial series. 

We should stress that the cross-overs corresponding to the radii of con- 
vergence of the Mayer and virial series do not necessarily correspond to sharp 
phase transitions, and may very well be determined by singularities off the 
positive axis. Nevertheless, they reflect changes in low-temperature physical 
behavior. This is formally analogous to resonances in quantum mechanics, 
when Green's function singularities off the real axis do not qualify as eigen- 
values, but can nonetheless affect the system's behavior. 



2 Setting 

We are interested in the statistical mechanics of a classical system of particles, 
in continuous configuration space, interacting via a pair potential v{\x — y\). 
Thus let V : [0, cx)) ^ M U {oo} and 

U{xi, . . . ,Xn) ■■= ^ v{\Xi-Xj\) 

l<i<j<N 

be the energy of a configuration of particles Xi, . . . , x^v G M'^. We assume 
that the energy is stable, i.e., there is a constant B > such that 

ViV G N, V(xi, . . . , XTv) e {R'^f : U{xi, ...,xn)> -BN. (2) 

In addition, we assume that the pair potential has finite range, i.e., v has 
compact support. For a given inverse temperature /3 > and A C M*^, the 
canonical partition function is 

Za(/3, N):=^ [ e-'^^^^^'-'^-Mxi ■ ■ ■ dxN, 

and the free energy per unit volume is 

/(/3,p) :=-lim^logZA(/3,iV). 

The limit is along oo, A = [0, with L ^ oo, N/L'^ p > 0. It is 

well-known that if the potential has no hard core (rhc = 0), the limit exists 
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and is finite for all p > 0; if the potential has a hard core, then for a suitable 
number p^p > (the close-packing density), the limit is finite for p < p^p, 
and infinite for p > p^p. Moreover the free energy p) is a convex function 
of the density p. 

The pressure at inverse temperature (3 and chemical potential ^ G M is 
p{/3,fi):= sup [pn- f{(3,p)). (3) 

0<p<pcp 

We call p(/3, p) the maximizer in the previous relation, if it is unique. Because 
of convexity, the density p(/3,p) is an increasing function of the chemical 
potential p. 

At fixed temperature, for sufficiently negative chemical potential, the 
pressure is an analytic function of the chemical potential and the activity z, 
with expansion 

Pp{(5,^) = z + Y,UP)z\ z = exp(/3p), (4) 

n>2 

the Mayer series, and the density is given by 

oo 

p(/3,p) = z + ^nhn{l3)z''. 

n=2 

Similarly, at low density, the free energy is strictly convex and analytic with 
expansion 

p) = p(log P - 1) - 5^ d^{P)p\ (5) 

n>2 

Eq. ([3]) gives, for p negative enough, 

Pp{^,p)=p-Y,{n~l)d^{P)p\ P = p{fi,p), (6) 

n>2 

the virial series. We would like to know how large z, or p, can be in those 
equations, and define 

R^^^((3) := supj^; > | Eq. Q is true with absolutely convergent series}, 
^"{(3) := sup{p > I Eq. (|5| is true with absolutely convergent series}. 

In principle, and R^"{P) can be smaller than the radius of conver- 

gence of the corresponding series: Eqs. (|4]) and (|5| might cease to hold before 
the series start to diverge. We do not know of any concrete example in our 
setting, but there is a simple type of situation where a similar phenomenon 
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arises in mean-field or Landau theories: suppose that a free energy F{m) is 
the convex envelope of some double-well potential, e.g., W{m) = (m^ — 1)^. 
Then F{m) = in |m| < 1, but at \m\ > 1 it becomes equal to W{m) > 0; 
in particular, F{m) ceases to equal its (trivial) Taylor expansion around 
before this expansion ceases to converge. 

For non-negative potentials, however, it is known [Ullll] that the domains 
of convergence coincide with the domain of equality of Eqs. (|4]), (|5]), and (|6]), 
so that in this case R^^'^{P) and R^"^{(3) are exactly equal to the radius of 
convergence. 

Furthermore we define 

Aisat(/3) := sup{/i e R I p{P, fi) is analytic in fi < ft}, (7) 
Psat(/3) := sup{i? G (0,pcp) I /(/3,p) is analytic in < p < i?}, (8) 

the chemical potential and density at the onset of condensation, i.e., the 
quantities associated with saturated gas. In the absence of a phase transition 
- for example, in one dimension -, Psat(/3) = Psat(/3) = oo. Another quantity 
of interest is 

p^^^(/3) := sup{p(/3,p) I exp(/3p) < (9) 

When there is no phase transition at = /J"^ log R^^^{(3), then p^^^{(3) = 

p(/3, p'^^^(/3)) . When there is a phase transition, the density may have a jump 
discontinuity and p(/3, is no longer well-defined; Eq. ^ states that 
in this case p^^^{(3) equals the left limit of p(/3,p) at p^'y(/3). 



i?^'' pMay p^^^ 







analyticity of f{l3,p) 




absolute convergence of J^^^ bnz{p)" 




absolute convergence of ^„ d^p" 







p 



Figure 1: Density axis, radii of convergence of the virial expansion (R™) and of 
the cluster expansion (p^'^^), and domain of analyticity of the free energy. 

We have the general bounds, illustrated in Figure [2| 

/3-Mog/2^^y(/3) < psat(/5), max(i?^'^(/3),p^^y(/3)) < psat(/3). 

With these notations, we can ask: 
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Question: When are the previous inequalities strict? When they 
are strict, is it nevertheless possible to give physical meaning to 
and T", even though in this case i?^"^ and do not 
correspond to phase transitions (i.e., points of non-analyticity) ? 

The main goal of this article is to show that the answer to the second question 
should be yes; moreover, in the presence of a phase transition for attractive 
potentials, the inequalites should be approximate equalities, in a sense spec- 
ified in the Corollary |3.3 and the Conjectures 3.4 and |3.9 below. 

We conclude this section with a description of the convergence criterion 
for the Mayer series that we shall use. Let 

En:= inf _U{xi, . . . ,Xn), = 0, (10) 
xi,...,xjve( 



be the ground state energy for N particles (without any volume constraint), 
and 

eoo := inf ^ = lim % < (11) 

(note that (En) is subadditive). The stability assumption on the pair po- 
tential ensures that Coo > —oo, and Eq. ([2| holds with B = — Cqo as optimal 
constant. We shall make repeated use of the following Theorem, which is a 
direct consequence of [121 Theorem 2.1], see also ^ for integrable potentials 
(without hard core). 

Theorem 2.1 (Mayer series estimates [IS]). Let v{\x — y\) be a stable pair 
interaction potential and rhc > the radius of the hard core. Set 



\\\v\\\ := |5(0,rhc)|+ / \vi\x\)\dx 

^IRd\B(0,rhc) 

where B{0,r^c) is the ball of radius centered at 0. Then 

R^'^'W) > -pTTTT^ (12) 

pe|||f III 

and for every < z < exp(/3eoo)/(/3e|||f |||), 

^n|6„(/3)|;z"-i < (e - 1)6"^'^°°. (13) 

n>2 

As an immediate consequence, we note that 

liminf /i,at(/3) > liminf /3-Mogi?^^y(/3) > e^. (14) 

Therefore every chemical potential /x < Cqo, as /3 — )■ oo, eventually falls into 
the gas phase. 
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3 Results and conjectures 

Our standard assumptions on the potential are the following: 

Assumption 1 (Minimal assumptions), f : [0, oo) — t- M U {oo} satisfies the 
following assumptions: 

• The energy is stable in the sense of Eq. (g. 

• f is everywhere finite except possibly for a hard core: there is a rhc > 
such that v{r) = cxd for r < The and v{r) < oo for r > The- 

• V has compact support, b := sup{r > | f (r) 7^ 0} < oo. 

• f is continuous in [rhcCxo), i.e., either f(rhc) = 00 and v{r) — )• 00 as 
r \ The, or t'(rhc) < c>o and f (r) — )• f (rhc) as r \ rhc. 

• V has an attractive tail: for suitable 6 > and all r G {b—6, b), v{r) < 0. 

Note that we allow for f (rhc) < 00, which is relevant for Radin's soft 
disk potential |T7j. The continuity of the potential is assumed in order to 
simplify statements on the low-temperature asymptotics. In particular, for 
sufficiently large volumes, 

lim hog Za{P,N) = -En. 

A similar statement holds of course without continuity assumptions, provided 
that we replace the infimum in the definition (10) by an essential infimum; 
the continuity allows us to avoid this measure-theoretic complication. 

Assumption [1] will be enough when working in the low density gas phase. 
For results that hold all the way up into a finite density region, we will make 
additional assumptions. We refer to every minimizer (xi, . . . ,Xn) G (M*^)^ 
of U{xi, . . . ,xn) as an A-particle ground state. Note that the attractive tail 
favors configurations where particles stick together. 

Assumption 2 (Ground state geometry and Holder continuity). For suitable 
a > 0, tq > 0, and every A G N, there is an A-particle ground state 
(xi, . . . , Xiv) e (M'^)^ such that 

• the interparticle distance is bounded below by tq: for alH 7^ j, \xi — 

Xj\ > ro; 

• the ground state fits into a cube of volume Na^: xi, . . . , xat G [0, N^^'^a]'^. 
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Moreover v{r) is uniformly Holder continuous in [ro, oo). 

The simplest example, in dimension two, of a potential satisfying As- 
sumptions [l] and [2] is Radin's soft disk potential [17] , which involves a hard 
core and a suitable attractive part. More general potential classes, again in 
dimension two, are given in f2T\. 

Assumption [2]is enough to ensure that various limits can be interchanged. 
In particular, if 

e(p) := lim 

p—>-cx) 

is the ground state energy per unit volume at density p, then 

. e(p) 
Coo = mm . 

o<p<pcp p 

Moreover e{p)/p has a minimizer p* < 1/a'^, i.e., the ground state has a 
finite preferred density. In [T71 [21] , p* is the density of particles in a simple 
hexagonal lattice. 

Our first result is about the low-temperature behavior of the Mayer co- 
efficients and should be contrasted with the alternating sign property for 
non-negative potentials [HI Chapter 4]. Recall the ground state energies En 
from Eq. (|To|. 

Theorem 3.1 (Mayer coefficients at low temperature). Suppose that v sat- 
isfies Assumption^ Then, for every fixed k, as (3 ^ 00, bk{(3) is eventually 
positive, and 

hm l3-'\oghW) = -Ek. (15) 

/3->c« 

Thus we may think of the Mayer series as 

/3p^5^z^exp(-/3E,). (16) 



The subsequent results are best understood with the approximate formula ( 16 ) 



in mind. We should stress that the approximation ( 16 ) can be derived with- 
out using Mayer expansions, see [6j for a quantum Coulomb systems result. 
Direct proofs are, in fact, much more instructive from a physical point of 



view; therefore Theorem 3.1 should be seen as a verification of the consis- 



tency of the Mayer series with the approximation (16). 

The next theorem builds upon a low temperature, low density result 
from [To] which we briefly recall. Suppose that v satisfles Assumptions [l] 
and [2] Then, for suitable /3o, po, Cq > and all (3 > f3o and p < po, 

Efe + rMogp 



<Copr'log/3. (17) 
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As explained in Appendix |B} po should be thought of as the preferred ground 
state density (an upper bound is Pq <l/a'^ with a as in Assumption [2|. The 
inverse temperature /3o is essentially determined by the condition exp(— < 
l/(a + RY, where 

z/* := inf(Efc - fccoo) > 0. (18) 
For potentials with an attractive tail, we have v* > Q |10j . 

Theorem 3.2 (Density increase around /i = Cqo)- Suppose that v satisfies 
Assumptions^ and^ and that x i— )■ v{\x\) is integrable in \x\ > r^c- Let 



Co, Po) Po > be such that Eq. (17) holds for all /5 > /3o and p < po. Then: 



• For every C > Co and suitable f3c > (3o > 0: 

V/3>/3c yp>e^ + Cp-HogP: p{p, p) > p,. 

• For every C > 1, all n & N and suitable I3{n, C): 

V/3>/3(n,C) Vp<eoo-C/3-Mog/3: p(/3, /i) < /?-". 

In particular, for every fixed /i > Cqo, as /9 — )• oo, the density is bounded 
away from zero, while for p < Cqo, it vanishes exponentially fast (Eq. (19) 



anticipates on Theorem 3.6) 



p> Coo ■ liminf p(/3,p) > po > 0. 

/3— >-oo 

p<eoo: p(/3,p) = 0(e-^"*). (19) 

Remark (Non-negative potentials). When f > 0, a similar change in the den- 
sity behavior occurs around p = 0, as the following two examples illustrate. 
For an ideal gas in continuum space, f3p = z, p = z, = 0. For a lattice 
gas with no interaction except the hard-core on-site repulsion, 

/3p(/3,yu) = log(l + 2;), p(/3,p) = — — , Coo = 0. 

As /3 — )■ oo, if p > is fixed, the density diverges (for the ideal gas) or ap- 
proaches the maximum density (for the lattice gas). For both the continuum 
and lattice gas, at fixed p < 0, the density goes to exponentially fast, but 



in contrast with the attractive potential case Eq. (19) there is no positive 
lower bound on the rate of exponential decay, u* = 0. 

A first consequence is an indication where the low temperature / low 
density solid-gas transition is located, if such a phase transition takes place. 
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Corollary 3.3 (Where to look for a solid-gas transition). Under the assump- 



tions of Theorem |3.2 , if Psat(/3) — )■ as /3 — )■ oo, then 

/Xsat(/3) = eoo + 0(r'log/3) 

as /3 — > cxD. 



The lower bound in the corollary follows from Eq. (12), noting that 



/^sat(/3) > /3 ^ log -R'^^y(/3). For the upper bound, suppose by contradiction 



that /isat(/3) — Coo ^ log/3 as /3 — t- oo. Then Theorem 3.2 tells us that 
Psat(/9) is bounded from below by some positive constant times po, contra- 
dicting the assumption Psat(/9) — ^ 0. Thus /isat(/9) — Coo = 0(/3~^ log/3). 

Remark (Lee- Yang theorem) . For a lattice gas on Z*^ with at most one particle 
per lattice site and attractive pair interactions v{x — y) < 0, the Lee- Yang 
theorem [HI Theorem 5.1.3] says that if there is a phase transition, then it 
must be at a chemical potential fi that satisfies 



1, 



i.e., /X = {^xf^Qv{x))/2. The right-hand side of the latter equality is readily 
identified with Cqo, the ground state energy per particle for the lattice gas. 
Thus Corollary |3.3 compares well with the Lee- Yang theorem. 



Remark (Widom-Rowlinson model). It is instructive to look at a contin- 
uum space model for which the existence of a phase transition is known, 
the Widom-Rowlinson model [22], see the review ^20j. Consider particles 
interacting via the energy 

Ua{xi,...,xn) = |AnU,^i5(x„l)| -iV|fi(0,l)| 

wher B{x,l) is the ball of radius 1 centered at x. The interaction is not a sum 
of pair interactions, but it qualifies nevertheless as an attractive, stable, finite- 
range interaction. The ground state energy per particle is e^o = — |-B(0, 1)|. 
An equivalent formulation is in terms of a two-species model with hard core 
repulsion between particles of different type: 



oo ^ 
Z 



E 



m 

N=0 



^2|A| 



e 

Ni,N2=0 



^ j,JVl iV2 /■ /■ 



10 



provided 

It is known [121 E] that for sufficiently high, equal activities zi = Z2, the 
system has a phase transition. In the one-species picture, a phase transition 
happens at low temperature and activity z = /3exp(— /3|i?(0, 1)|), or chemical 
potential 

/i = -|5(0, 1)1 + /3-Mog/3 = Coo + /3-Mog/3. 



Again, this matches Corollary |3.3 



A second consequence of Theorem 3.2 is that, even when there is no phase 



transition - for example, in one dimension -, there is nevertheless a change in 
physical behavior around ^ = Coo'- consider the family of curves /i p{l3,fi) 
around /i = Coo- At /3 = oo, it has a jump of size > po- At (3 large but 
finite, there could be either a jump, or the curves resemble the occupation 
numbers of fermions around the Fermi energy. Hence there is either a phase 
transition, or a fast increase from small to large density. 

We would like to propose this as a a physical interpretation to the domain 
of convergence of the Mayer series, for attractive potentials, based on the 
following conjecture: 

Conjecture 3.4 (Mayer series' radius of convergence). Suppose that v sat- 
isfies Assumptions [T] and |2j Then 

lim rMogi?^^y(/3) = eoo. (20) 
Note that for pair potentials whose finite part is integrable, we have the 



lower bound Eq. ( 14 ) on the liminf . Hence the only part that is open in the 
previous conjecture is an upper bound on the limsup. 

In fact, a rigorous statement is available for the radius of convergence of 



the finite volume Mayer series. It is proven by combining Theorem |3.1| with 
the bounds from [H]. First recall that the pressure (3p\{(3, z), defined via the 
logarithm of the grand-canonical partition function in a finite box A = [0, L^, 
has a Mayer expansion similar to Eq. Q, with volume-dependent radius of 
convergence R^"^{(3). Note that liminf ia^oo < see [H 

Eq. (4.2)]. 

Theorem 3.5. Let the pair interaction satisfy Assumptions\l\ and\^ Then, 
if we let first |A| — )■ oo along cubes, and then /3 — )■ oo, 

lim limsup/3~^ logi?^^^(/3) = lim liminf log -R^^^(/3) = Coo- 

/3^oo |A|^oo |A|-^oo 
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We are now heading towards similar interpretations for the virial expan- 
sion. First, however, we need to understand better the gas phase /i < Cqo- 
The reason is that inside the gas phase, there might be "chemical" transi- 
tions [9j, for example, from monatomic to diatomic gas. When this happens, 
the radius of convergence of the virial series is determined by that cross-over, 
and the virial series ceases to converge before any sharp phase transition is 



observed, see Theorem 3.8 



The next theorem is a grand-canonical version of results from [TU] and 
should be compared to the atomic or molecular limit for quantum Coulomb 
systems [31 EJ [7j. (See also a result for the classical one-dimensional two- 



component plasma [12] ■) Recall the quantity z/* > defined in Eq. (18) 



Theorem 3.6 (Possible cross-overs inside the gas phase). Suppose that v 
satisfies Assumption^ and that v is integrahle in \x\ > rhc- Then for every 
fixed < Coo, 

lim (3-Hogp{l3,fx) = -inf(Efc - kfx) < -u* < 0. (21) 

If in addition {Ek — kfi)k£m has a unique minimizer k{fj,) G N, then as (3 ^ oo, 

M/3,/i) = ^^(1 + 0(1)). (22) 

The interpretation is that the gas is, approximately, an ideal gas of 
molecules consisting of /c(/i) particles each, with effective activity z'' exp{—(3Ek), 



see also Eq. (16). An illustration with coexistence curves in the density- 



temperature plane is given in [TUl Figure 4]. 



Remark. By now we have two auxiliary variational problems: in Eq. (21) 



we minimize — kjj, with respect to k at fixed /z; in Eq. (17), we minimize 
{Ek + (3~^ log p)/k with respect to k at fixed /3 and p. In Appendix [A| we 
show that these two minimization problems are equivalent and discuss their 
properties. 



Theorem 3.6 covers two different scenarios, depending on the value of 



II, := inf (23) 

k>2k -1 

The quantity pi separates a region dominated by monomers from a region 
dominated by larger groups of particles: for fixed p < pi, Ek — kp has the 
unique minimizer k{p) = 1, and for p > pi, every minimizer is > 2, see 



Lemma A. 3 As a consequence, for sufficiently negative chemical potentials 



we observe a monatomic gas {k{p) = 1). If /xi = Coo, this is all we see in the 
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gas phase. If /Ui < Coo, as we increase the chemical potential, we observe a 
transition from monatomic to polyatomic gas before the gas condenses. 

The existence of such a transition becomes very natural when we look 
at a concrete example, taken from [5, Sect. 6]. Consider a pair potential 
with a hard core and two potential wells, a deep well at small distances, 
and a shallow well at larger distances, separated by a repulsive {v > 0) part 
at intermediate distances. The deep well favors small groups of particles 
(pairs, triangles or tetraeders, depending on the dimension), arranged at 
larger distances determined by the shallow well. We may think of a solid 
made of molecules instead of atoms. It is natural, then, that the solid forms 
after atoms gather in molecules. A rigorous statement with a proof of ni < 
Coo, in dimension one, can be found in [5]. 

The previous example suggests a relationship between the geometry of 
ground states and the existence of a cross-over inside the gas phase. An 
interesting open question is, therefore, whether the conditions from [211 C] 
ensuring a crystalline ground state with hexagonal lattice (one particle per 
unit cell) also imply /ii = Coo- A much weaker result is the following: 

Proposition 3.7 (Sufficient criterion for the absence of polyatomic gas). Let 
V be a stable pair interaction with attractive tail. 

1. If for all m,n ^N, 

Em+n+l < Em+1 + En+1 (24) 

then /ii = Coo- 

2. Suppose that v{r) has a hard core rtc > and v{r) < for r > rhc- 



Then, in dimension d = 1, the inequality (24) holds for all m,n, and 
we have fii = e^o- 



Eq. (24) should be read with a "gluing" operation in mind: instead of 
juxtaposing m and n-particle configurations in space, as is done in order to 
derive the subadditivity Em+n < E^ + En, we glue two configurations with 
m + 1 and n + 1 particles in one point. 

Remark. We owe to G. Friesecke the following remark: in statement 2. of 



Prop. 3.7, we may replace the assumption that v has a hard core by a state- 
ment of the type "f is sufficiently repulsive near the origin" , formulated for 
example through inequalities for derivatives of the potential. Precise state- 
ments are, already in dimension 1, surprisingly involved. 

After this excursion into ground states, let us come back to the virial 
expansion and compare the density of saturated gas psat of Eq. (Is]) with the 
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virial radius of convergence R"'" and with p^^^ defined in Eq. g. Recall 
u* > from Eq. (18) and let Ui := — /ii with /ii as in Eq. (23). The quantity 
Ui is a canonical version of the grand-canonical threshold fii. We note that in 
general z/i > i/*, and ui > v* if and only if /ii < Coo, i-e., if and only if there 
is a monatomic-polyatomic cross-over inside the gas phase (see Lemma A.2). 

Theorem 3.8 (Comparison of i?^"", p'^^^, psat and vi). Suppose that v satisfies 
Assumptions\^ and^ Then 



liminf /3-Mogpsat(/3) > liminf /3-Mogp^'y(/3) > - 



/3— >oo 



liminf /3-Mogi?^'^(/3) > -Vi. 

If in addition pi < Cqo and Ei^/{k — 1) has a unique minimizer, then 

-z/i < 



lim /3-Mogi?^'^(/3) 
and as (3 -^oo, R"'{P) < p^^y(/3) < psat(/3). 



-z/ 



(25) 
(26) 

(27) 



Eq. (27) tells us that if there is a monatomic-polyatomic cross-over inside 
the gas phase, then that cross-over determines the radius of convergence of 
the virial expansion, and the virial expansion ceases to converge before there 
is any phase transition. 

This result should hold without the technical assumption that Ek/{k — 1) 
has a unique minimizer. In fact, it is natural to think that Eq. (27) extends 
to the case ui = v* , so that the radius of convergence of the virial expansion, 
for attractive potentials, is always determined by the first cross-over - either 
from monatomic to polyatomic gas, or directly from small density, monatomic 
gas, to large density; the latter cross-over possibly being a phase transition 
(in d > 2). We also have a conjecture on the behavior of p^^^ and psat 



analogous to Conjecture 3.4 and Corollary 3.3 



Conjecture 3.9. For interactions with an attractive tail. 



lim /3-Mogp^^y(/3) 

/3— >co 



If in addition there is a low-density, low-temperature phase transition, i.e., if 
Psat(/3) — )■ as /3 -7- oo, then 

lim p-'\ogpUP) = -v\ 

/3— >co 
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Let us recall that the line p = exp(— /3z/*) has the following physical 
interpretation, proven in [10]: at densities that are very small but higher 
than exp(— /3z/*), particles tend to gather in very large clusters (i.e., groups 
of particles close in space), even though the system is dilute. At densities 
smaller than exp(— /3z/*), particles stay for themselves or form small groups 
- this is the gas phase discussed above. 

Finally, we have partial results on the low-temperature asymptotics of 
the virial coefficients dk((3) from Eq. ([s]), to be compared with Theorem 3.1 



en 



Proposition 3.10 (Virial coefficients in the absence of polyatomic gas). Let 
V satisfy Assumption^ Suppose that Eq. (24) holds for all m,n eN. Th 
/^i = Coo and for all k > 2, 

limsup/3"Mog4(/3) < -Ek. 



If in addtion the inequality (24) is strict for all m,n E N, the previous 



inequality for the limsup becomes an equality for the limit. 

Proposition 3.11 (Virial coefficients in the presence of a monatomic-di- 
atomic transition). Suppose that v satisfies Assumption^ and in addition 
Pi < God and Ek/{k — 1) has the unique minimizer k = 2. Then for every 
k > 2, as P ^ oo, dk{P) eventually has the sign of (—1)^, and 

lim r'log((-l)'4(/3)) = -ik-l)E2 > -Ek. 

In particular, \dk{P)\ / hk{l3) — )■ oo as {3 ^ oo. 



For the two-well example from p. 13, we expect pi < e^o and Ej^/{k — 1) 



should have the unique minimizer p = ci + 1, with d the dimension of the 
configuration space. In one dimension, this is proven [5l Sect. 4], and gives 



an example to which Prop. 3.11 applies. For higher dimensions, we note that 



the natural generalization of Prop. 3.11 when pi = Ep/{p — 1) for a unique 
p > 3 is 

rfi+r(p--i)+g(/3) = i-iydq+i{l3)exp(-l3r{Ep + o{l)), 

r E No, q = 0, 1, ... ,p — 2. We leave the proof, or disproof, as an open 
problem, and do not exclude surprises - it is not impossible that additional 



conditions, in the spirit of Eq. (24), are needed. 

In any case, we see that the behavior of the virial coefficients is more 
complex than that of the Mayer coefficients. In the absence of a cross- 
over, at low temperature, each virial coefficients is eventually positive - as 
envisioned by Mayer and Mayer [13], Chapter 14(f)]. Note, however, that 
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in [13] the authors argue that there is a temperature below which all virial 



coefficients become positive; this statement is much stronger than Prop. 3.10 



When there is a cross-over, in contrast, the sign of the coefficients varies in 
a periodic way. 



4 Mayer coefficients at low temperature 



Here we prove Theorem 3.1 We use the usual short-hand Vij = v{\xi — Xj\), 
and fij as in 

exp{-l3v{\xi - Xj\)) = exp{-l3vij) = 1 + fij. 
We recall the expression of the Mayer coefficient: it is known that 

W) = hy2 [ T] f^Jix)dx2■■■dxk, xi:=0. (28) 



7 conn. v» I (ij)&'y 

The sum is over connected, undirected graphs 7 = (V, E) with vertices 
l,...,/c, and 11(4^)67 product over edges E E, i < j (no self- 

edges (a))- 

Let us start with a look at the /3 — ?■ 00 behavior for an individual graph. 
Observing that 



we get 



(1 + 0(1)) exp{-f3vij{x)), Vij{x) < 0, 
-1 + 0(1), Vij{x)>0, 



n |i^,(a.)| = (l + o(l))exp(-/3 J2 %N1K(^)<0)). 

(«i)G7 («i)e7 

In the exponent, only negative interactions appear. As a result, we may end 
up with energies much smaller than the ground state energy, seemingly con- 



tradicting Theorem |3.1 1 The reason is, of course, that there are cancellations 
between different graphs. In order to get a hold on them, it is convenient 
to do separate book-keepings for "positive" and "negative" edges. Given 
X = (xi, . . . , Xk), we define 

£+{x) := {{i,j} \1 <i < j < k, Vij{x) > O} 
S-{x) := {{i,j} \ l<t<j <k, Vij{x) < 0}. 

and let 'j^{x) be the graphs with vertices 1, . . . ,k and edge sets S^{x). 
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The next simplifying observation is that if the interaction has a finite 
range R > 0, fij{x) vanishes as soon as |xj — Xj\ > R. Therefore we define, 
for X = (xi,. . .,Xk), 

£{x) := {{'i, j} \ l <i < j <k, \xi- x,j\ < R], 

and let '~i{x) be the graph with vertices 1, . . . , and edge set £{x). We call a 
configuration x connected if the graph "^{x) is connected, and write lconn(a;) 
for the corresponding characteristic function. With these notations, for every 
configuration x and every graph 7, 

n k{^) 7^ ^ £(7) c 8{x), 

(«i)e7 

and if 7 is connected, so is cc. 

We are going to compare the Mayer coefficient with a partition function 
for connected configurations, 

Zt{P) ■■=T, f e-''^(°'^^'-'^'=)leonn(0, X^, • • • , X,)dX2 " " " dx,. (29) 

J(Rd)fc-l 

Lemma 4.1 (Cluster partition function vs. Mayer coefficient). 

- / Y\ fij{x)lconn{x)dX2- ■ -dXk, Xi = (30) 



7 not conn. (^ij^^^ 



where the sum extends over graphs 7 with vertices {1, . . . , /c} that are not 
connected. 

Proof. In the integral for Z^^(/3), write as usual exp(— /3t>jj) = 1 + fij and 
expand. This gives a sum over graphs. The graphs that are not connected 



correspond to the right-hand side of Eq. (30). The connected graphs yield an 



integral similar to Eq. (28), except that there is the additional characteristic 
function Iconnix). Noting that Wi^ij)^^ fij{x) vanishes if 7 is connected and 
X is not connected, we can drop the characteristic function without changing 
the value of the integral, and obtain Eq. ([30|). □ 



For 7 a graph with vertex set {1, . . . , /c}, and x = {xi, . . . , Xk) G (M'^)'^ a 
configuration, write ^~{x) (1 7 for the graph with vertices 1, . . . ,k and edge 
set £(7) n £~{x). Thus 7^(cc) fl 7 is the subgraph of 7 consisting of the 
negative edges. 
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Lemma 4.2. Let k eN and 7 a graph with vertices 1, . . . , /c with connected 
components of size ki, . . . ,kr, r & N, ki — k. Then 



E n 

7: 7 (a;)n7=7 (ij)67 



for some suitable Ck > which does not depend on /3 or r, 

similar estimate holds, for r > 2, if the sum is further restricted to graphs 7 

with 7~(a;) 07 = 7" that are not connected. 

Remark. The lemma is also true for a sum further restricted to graphs that 
are connected. It becomes wrong, in general, for doubly connected graphs. 

Proof. Consider first the case r = 1, i.e., 7" connected. Then 



E n z^^^^) 

7: 7-(x)n7=7- (ij)ej 




E n z^^^^) 

ece+{x) {ij)ee 



n 



-pVij(x) 



Noting that for a negative edge, < fij < exp(— /Jvjj), it follows that 

0^ E W h{x)<exv{-p[ E %(^)+ E ^^^H) 

7: 7-(a;)n7=7- (ij)e7 (u)e7+(a:) 

< exp(-^C/(xi,...,Xfe)) < exp(-;5£;fc). 

Next, suppose that r > 2 and that we restrict the sum to graphs 7 that 
have 7~(cc) 07 = 7" and are disconnected. Such a graph can be constructed 
from 7~ in two steps: first, add positive intra-component edges, i.e., edges 
{ij) that have vij > and connect two labels i.j belonging to the same 
connected component of 7". There is no restriction on the choices of such 
edges. Second, add positive inter-component edges. There is a restriction 
on how many edges we may add, since the resulting graph is required to be 
disconnected. Let T be the collection of allowed inter-component edge sets. 
The sum to be estimated takes the form 



E n n 



( n /^.•(-)) ( n ^"'"^^^0 



. {y)e7g 



{ij)e7q (as) 



(31) 
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Here , . . . , 7^ are the connected components of 7 , and ^^{x), . . . , 7^(ic) 
have as edge set the positive intra-component edges. Noting that — 1 < fij < 



for every positive edge, we deduce that (31) has absolute value bounded by 



1^1 exp(-/3(Efc, + + 

Since \J^\ can be bounded by some /c-dependent constant, independent of 7", 
this concludes the proof. □ 



Proof of Theorem 3. 1 Let x = {xi, . . . ,Xk) be an arbitrary configuration 
and 7 a graph with vertices 1, . . . ,k that is not connected. Then '~f^{x) fl 7 
is not connected either. Therefore 



7 not conn. (ij)Siy 



E En/. 

7~ not conn. 7 not conn.: (ij)^'y 

7~ {a;)n7=7~ 



Lemma 4.2 then yields a bound on the absolute value of the form 
Cfcexp(-/3(Efc, + ■ ■ ■ + EfeJ) < Cfeexp(-/3(r - l)e) exp(-/3Efc) 



(32) 



Here we have used that for potentials with an attractive tail, for suitable 
£ > and all k,q G N, -Efe+g < + ~ e (see the appendix in [TU]). 
Since the set of connected configurations (0, 0:2, ... , x^) has a finite Lebesgue 



volume, the integral on the right-hand side of Eq. (30) has an upper bound 



similar to the right-hand side of Eq. (32). On the other hand, because of 



the attractive tail of the potential, ground states are always connected. The 
continuity of the potential therefore yields 

lim /3-MogZ^'(/3) = -Efc, 



and we conclude from Lemma 4.1| that for every fixed k, as /3 — )• 00, 

6,(/3) = (1 + 0{e-'^))Ztm = exp(-/3(Efc + o(l))) . 



□ 



5 Virial coefficients and absence of polyatomic 

gas 

The virial series coefficients (i„(/3) have an expression as a sum of integrals. 



indexed by graphs, similar to Eq. (28); the sole difference is that the sum is 



over doubly connected graphs only (see, e.g., [9]). Thus we may write 
bn{P) = dn{/3) + a sum over graphs that are not doubly connected. 
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By Lemma 4.1, the main contribution to bn{P) comes from integrals around 



the ground state. In dimension 2, the ground state could resemble, for exam- 
ple, a hexagonal lattice and should be doubly connected; thus one might think 
that the main contribution to comes from doubly connected graphs and 
therefore = (1 + o(l))f/n(/3). 

It turns out that this naive guess leads to the right answer under the 



assumptions of Prop. 3.10 , but fails in the setting of Prop. 3.11 In the latter 
case, both and the sum over graphs that are not doubly connected are 

much larger, in absolute value, than 6„(/3). 

The reason for this complex behavior is that the cancellations between 
different graphs become rather subtle to handle; see also the remark after 
Lemma |4.2[ Therefore we work instead with an expression of the virial 



coefficients in terms of Mayer coefficients. It is known [HI Chapter 5, Eq. 
(25.30)] that 



dnW) 



E' (-1) 



(Earn,)-! 



(33) 



m2,. 



with a sum over (m2, . . . , m„) G Nq such that 



and 



a(m2. 



^(j - l)mj = n - 1, 
i=2 

^ _ (n-2 + ^2"^j)' -QJ 



> 0. 



j=2 3 



(34) 



(35) 



One may check that the contribution of the vector with m„ = 1 and m2 
■ ■ ■ = m„_i = is equal to bn{(3). 



The form of Eq. (33 ) becomes very natural when the virial series is derived 



directly with the help of a cluster expansion in the canonical ensemble, as 
recently done in [16], see Appendix [Cj 



Lemma 5.1. Suppose Eq. (24) holds for all m,n G N. Then, for all 
(ma, . . . , m„) G N^J"^ satisfying Eq. (|34|), 



i=2 



(36) 



If the inequality (24) is strict for all m,n, then the previous inequality is 
strict too. 
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Proof. We proceed by induction over r = Yyl^j- -^^^ r = 1, the inequality 
is a trivial equality. For r = 2, the inequality (36) is the same as (24), and 
true by assumption. For the induction step, suppose that the statement is 
true, at r, for all n G N. Let (m2, . . . ,mn) satisfy Eq. (34) and such that 
YI2 "^j = + 1 > 2. Write mj = m'j + Sjk for some A; G {2, . . . , n}. Then 



En = Ei+(j22{j-i)m'.)+{k-i) < -^i+E2(i-i)K + -^'^ ¥S 
Ek 



i=2 



i=2 



which proves the claim. The procedure for strict inequalities is exactly the 
same. □ 



Prop. 3.10 is an immediate consequence of Eq. (33), Lemma 5.1 and 



Theorem 13.1 



Proof of Prop. (3.10). By Lemma |5.1 and Theorem 3.1, all terms in the 
sum (33) are of order at most exp(— /3i?„(l + o(l)). If the inequality (24) is 
strict, the dominant contribution comes from (r = 1), which is equal 

to exp(— /3i?„(l + 0(1))), again by Theorem 3.1 □ 



Lemma 5.2. Suppose that fii = E2 < e^o and E^/ {k — 1) > E2 for all A; > 3. 
Then, for all (m2, . . . , m^) G Ng"^ satisfying Eq. ([34|), 



[n 



1)E2<J2^,E, 

i=2 



with equality if and only if m2 = n — 1, 



nir. 



0. 



Thus the main contribution to hn{P) — dn{fi) comes from graphs whose 
doubly connected components all have size 2. 

Proof. By assumption, Ek > {k — 1)E2 for all A; > 3, from which we obtain 

r 

Ek, + --- + Ek^> Y^iki - 1)E2 = {n- 1)E2 
1 

for all ki, . . . ,kr > 2, with equality if and only if all fcj's are equal to 2. 
Writing rrij for the number of i's such that ki = j, we obtain the desired 
statement. □ 



Proof of Prop. 3.11 By Lemma 5.2 and Theorem 3.1, in Eq. (33) all terms 
are negligible except the one for m2 = n — 1, = ■ ■ ■ = nin = 0. □ 
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We conclude with the proof of the sufficient criterion for the absence of 
polyatomic gas. 



Proof of Prop. 3/7. 1. Eq. (24) says that n ^ En+i is subadditive, thus 



/ii = inf 



E. 



n+l 



lim 



n 



2. Suppose d = 1 and v{r) < for all r > The and v{r) = oo for r < rhc- 
Let m,n eN and x, y he m + 1 and n + 1-particle ground states. Without 
loss of generality we may assume that a;i < ■ ■ ■ < Xm+i, Z/i < ■ ■ ■ < Vn+i, and 
2/1 = Xm+i- We construct a m+n+l-configuration z by gluing the two ground 
states: set Zi := Xi,..., Zm+i '■= Xm+i and -Zm+j 1/j for j = 1, ...,« + 1. 
Then 



m+n.+l 



<f/(2) 



E, 



m+l 



+ + (37) 

j=i k=i 



Because of the hard core, all particles have mutual distance > rhc and neg- 



ative interactions. Hence the double sum in Eq. (37) is smaller or equal to 
zero, and we deduce Em+n+i < E^+i + En+i- Since m and n were arbitrary, 
applying the sufficient criterion from 1., we get pi = e^o- □ 



6 Bounds for the density p(/5, /i) 



In this section we prove Theorems 3.2 and 3.6 



Proof of Theorem 3.2. 1. Suppose > Coo + 0(3 ^ log/3 with C > Cq. Write 



p = exp(— /3z/) and p{iy) = inffcgN(-Efc — i')/k. For v < u* , p{u) = Coo, and for 
u > u*, p{ij) < Coo, see Appendix [A} Using Eq. (17), 

p(/3,/i) = sup {pp- f{(3,p)) 

0<p<pcp 

> sup {pp-f{(3,p)) 

0<p<po 



> sup ipp-~p{iy)p — Cop(3 ^ log l3 

0<p<po ^ 



> sup ((/U - Coo - Co/3 Mog/3)p 
= po{p - Coo - Co/3"^ log 13) , 
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since > Coo + Co(3^^ log (3. On the other hand, let p = p{f3,fi) be any 
maximizer of pfi — f{f3,p). If p > po, we are done. If exp(— /3z/*) < p < po, 
then 

p{p - Coo + log /3) > p) > Po (p - Coo - C'o^^^ log /3) , 

which gives 

^<e^ + Co^^r'log/3. 
Po- P 

Since p > eoo + C/?""*^ log/3, we obtain Co(po + p) > C(Po ~ p) whence 

^ C-Cq 

Thus we are left with the case p < exp(— i.e., > v* . Noting p{v) > 
Coo — for all z/, we get 

pp - p(eoo - z^) + CoP/3"^ log (3 > p{l3, p) > po (p - - Co/3~^ log /3) . 

Since pz/ = i^e"^"^ < we obtain 

r' + c^o(po + p)r'iog/3 



P < Coo + 



< eoo + 



Po -P 

r' + a(po + e-^'^*)r'log/3 



Po - e- 



eoo + Co(l + o(l))r' log/3, 



which for sufficiently large /3 is in contradiction with the assumption on p. 
2. For p < Coo — C/3~^log/3 with C > 1 we use Theorem 3A Define 
> by 

/2e-^^-/3|||t;||| = 1/e, (38) 



see Eq. (12). For K gN, using Eq. (13), we have 

K 



p(/3,p)<5^A;exp(/3(A;p-E, + o(l))) + (|) ;^ ^ A;|6,(/3)|/?^-^ 



fc=i 

<- g-/3(i.*+o(l)) 



k=K+l 



R 



(e — l)ze 



+ (e - 1) exp (^ir [(-C + l)/3-^ log /3 + /3-i(l - log 1 1 1 1)] ) 
As /3 — )■ oo, the second term is of order Since K could be 



chosen arbitrarily large, this completes the proof of Theorem 3.2 



□ 
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Proof of Theorem 3j6_. We proceed analogously to the proof of the second 
part of Theorem 3.2 Let R = exp(/3(eoo + o(l))) be as in Eq. (38). Fix 
fi < Coo- For K eN, using Eq. ( [I3| ), 

K 



- ^ A;exp(/3(A;^ - + o(l))) 



k=l 



< 



R 



K 



l)ze- 



(39) 



Since Ek — kfi = k{eoo — + o(l)) — )■ oo as fc ^ oo, there is a finite k{fi) 
minimizing Ek — kfx. Choosing K large enough so that K > k{fi) and K{eoo — 
fi) > I mik{Ek — kfi)\, Eq. (21) follows from the inequality (39). 

If {Ek — A;/i)fceN has a unique minimizer G N, the previous argument 
actually yields 

p(/3,/x) = fc(/x)6fc(^)(/3)e'^^('^)'^(l + o(l)). 
An analogous argument gives 



and Eq. (21) follows. 



□ 



7 Asymptotics of R^^^, p^at, P^'^ and R"'' 



Here we prove Theorems 3.5 and 3.8 



Proof of Theorem \3.5\ If the activity satisfies Eq. (12), then z < R^^^{P), 
see [151 Theorem 2.1], and the lower bounds follow just as in Eq. (14). For 
the upper bound, we use a result from [Tl]: for all A; G N, 



/ A:exp(-/3eoo) yA^-D 



Here 6fc,A(/3) are the coefficients of the finite volume pressure-density series. 
They converge to bk{f3) as |A| — > oo, whence 

r pMay.^N ^ /^ /cexp(-/3eoo) y/(fe-i) 

hmsupi?^ ^ /3 < (— — 1 
lAKoo V(fc - 1) 



We deduce from Theorem 3J^ that for every k eN, 

limsuplimsup/3~^ logi?^^^(/3) < — — ^ 

/3^>oo |A|->-oo k — 1 



We conclude by letting /c — )■ cxd in the upper bound. 



□ 
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For the proof of Tfieorem 3^, we start with the lower bound on the den- 
sity of saturated gas Psat(/3) and the density p'^^^(/3) dehmiting the physical 
parameter region covered by the Mayer series. 



Proof ofEq. (|25|. We observe that psat(/3) > P^'^(/3) > p(/3,p) for all /3 and 
all p < Psat(/3)- By Eq. (14) and Theorem 3.6, it follows that for every 

P ^ Coo) 

liminf/3-Mogpsat(/3) > liminf /3-Mogp^^y(/3) > sup(A;p-Efc). 



Noting that 



sup sup(A;p — Ek) = sup sup (A;p — E^) = sup(A;ec 

^J.<eao fcGN fceN ^<eoo fcSN 



we deduce Eq. (25). 



□ 



Proof of Eq. ( 26 ) . Write the pressure-density series as Pp = p+X]n>2 <^n{/3)p^' 
We start from the contour integral, see |TT] or |T8l Chapter 4.3], 



c„(/3) 



dz 



27ri Jq nzp{z) 



n-i ■ 



Here the density p{z) = Xl^i is extended to complex activities 
z, and we integrate on a circle C of radius exp(/3p) < exp(/3pi). For /3 
sufficiently large, we know that exp(/3p) < R^^^{P), and we are going to 
check that \p{z)\ > for \z\ = exp(/3p). To this aim we write 



k>2 



For every fixed K > 2, as P ^ oo, 

K K 



J^klhWW)'-' = 5^fcexp(/3(A;- l)[p- ^ + o(l)] 

k=2 k=2 

K 

< Y,ke^v[P{k -l){^-p, + 0(1))) 

k=2 

< const(-ft', p) exp^/3(p — pi + o(l))) . 

Furthermore, if = exp(/3(eoo + o(l))) is as in Eq. (38), and (3 sufficiently 
large so that exp(/3p)/i? < exp{—f3e) with suitable e > 0, 



< 



R 



K 



(e - l)exp(-/3eoo) 



< (e-l)exp(-/3(ir£ + eoo)). 
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We choose K G N such that Ke + Coo > /^i — /U and combine the previous 
estimates. We obtain that as /3 — )■ oo, J2k>2 ^\^k{.P) \ \z\'^~^ is of order at most 
exp(/3(/i — yUi)) and, in particular, goes to 0, so that p{z) 7^ 0. 

We can plug the lower bound for \p{z)\ into the contour integral. This 
yields 

1 



\Cn\ < 



n 



[(l + o(l))exp(/3/x 
whence > (1 + o(l)) exp(/3/i) and 



liminf/3-Mogi?^'^(/3) > /i. 



n-1 



This is true for every fi < pi = —ui, and the inequality (26) follows. 



□ 



Proposition 7.1. Under the assumptions of Theorem \3.S\ suppose that in 
addition fii = E2 < Coo- Then, for suitable (3i,Ci,6 > and all P > the 
equation 



k=2 



has a solution Zo{f3) at distance < Ciexp{—ef3)/b2{f3) 0/ —1/(462 o,nd 
p{zo{(3)) 7^ 0. // pi = Ep/{p — 1) < Coo for some p > 3, a similar state- 
ment holds with — 1/462 replaced with one of the roots of the equation 
l+p%p{l3)zP = 0. 



Proof of Eq. (27). For sufficiently small p, the density-activity relation can 
be inverted: there is a function C{p), analytic in a domain containing 0, such 
that for small z, ({p{z)) = z, and the restriction of ( to some neighborhood 
of is injective. The virial series is given by the composition 



00 00 



n=l 



n=l 



For sufficiently small z, we have ({p{z)) = z and 



P'{p{z)) 



zp'{z)' 



(40) 



The relation extends by analyticity to every domain D such that p is analytic 
in 2 G -D and P is analytic in p{D). Now, from Prop. 7.1, we know that 
p'{zq)/ p{zq) = with < l^ol < R^^^- Eq. cannot be true at 2; = zq- 
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Let D be an open disk centered at with radius > |zo|- The function P 
cannot be analytic in all of p{D), hence there must be some z E D such that 



R''"<\piz)\<J2k\bkm^ 

k=l 



(41) 



Let S > small enough so that fii + 6 < e^o- For sufficiently large (3, we may 
choose \z\ < exp(/3(/ii + 6)) in Eq. (41). The usual procedure shows that 

limsup/3-Mog/2"'^(/3) < sup(A;(/ii + 5) - E^) . 

Noting that sup f,{kfi — E^) is locally bounded and convex, hence continuous, 
in /i < Coo, we can let 5 \ 0, which yields 



limsup/3"Mog/?^'^(/3) < snp{kpi~Ek] 



-Z/1. 



□ 



Proof of Prop. \ 7.1\ The idea is to use an implicit function theorem, perturb- 
ing around exp(— /3) = 0. We give the proof for pi = E2 < 00. The proof for 
/ii = Ep/{p — 1), p > 3, is similar. For sufficiently large /3, &2(/3) > 0. Let 



ak{(3) := k 



2 W) 



z := h2{P)z. 



The equation to be solved becomes 

l + 4i + ^afc(/3)z'=-^ = 0. 

fc>3 

By assumption, for suitable A > and all k > 3, Ek > fii + {k — 1)A, and 
Hi = E2 < Coo- Therefore, for every fixed > 3, as /3 — )■ cxd, 

|afc(/3)| < A;2exp(-/3(A;-l)(A + o(l))) ^0. 

Moreover, for every s > and K G N, 



k=K+l 



1+s 



h2{PY-^ 



< ( sup(A; + l)^(62(/3)i?)-M(e - 1)6"''"° 



with R = exp(/3(eoo + o(l))) as in Eq. (38). Thus 

b2iP)R = exp(/3(eoo -E2 + o(l))) ^ 00. 
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It follows that for every s > 0, and suitable > 0, as /3 — t- oo, 

oo 

=0(exp(-£,/3)). 

fc=3 

Now let X be the Banach space of sequences {ak)k>3 with weighted norm 
I |a| I := J2k ^^Wk\- Set 

F{a, z) -.= 1 + 42 + ^ akZ^~^. 

k>3 

The equation F{0,z) = has the unique solution zq = —1/4, and in a 
neighborhood of (0,^0)5 for suitable choice of s, F is continuously Frechet- 
differentiable. Moreover ^^^(O, z) = A ^ 0. As a consequence, we can 
apply a Banach space implicit function theorem |23[ Chap. 4]. It follows 
in particular that as ||a|| — 0, the equation F{a,z) = has a solution 
z{a) = ^0 + Odlall). Applying this to a{(3) = {ak{f3))keN, we obtain the 
solution b2il3)zoi(3) = -1/4 + 0{e->^'^). 
For the density, we observe that 

pizom = zoif3) (1 + 262(/3)^o(/3) + 0(6"^^^)) = -^^o(/3)(l + 0(e-^^^)). 
It follows that for sufficiently large (3, p{zo{f3)) 0. □ 



A Two auxiliary variational problems 

Throughout this section we assume that w is a stable pair potential with 
attractive tail. Consider the following two variational problems 

z/(/i) := inf {Ek - kp), fi < e^c, 
fceN 

ii(h') := inf — ^- , u > 0. 

km k 



The first variational problem appears in Theorem 3.6, and minimizers k{p) 



correspond to the favored size of molecules in the gas phase as /3 — )■ cx3 at 



fixed /i. The second problem appears in Eq. (17) and, as shown in [10] . 
minimizers k[i') correspond to favored cluster or molecule sizes as /3 — )■ 00 
and p — )■ along p = exp(— at fixed z/. 

Recall that for potentials with an attractive tail, u* := mik{Ek — keoo) > 0. 

Lemma A.l (Concavity, monotonicity and equivalence). Let f be a stable 
pair interaction with attractive tail. Then: 
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1. The function fi h-)- z/(/i) is strictly decreasing, piecewise affine and con- 
cave in G (—00,600]- The function u (-)■ /i(z/) is decreasing, piecewise 
affine and concave in z/ G [0, 00); it is strictly decreasing in z/ G [z/*, 00) 
and equals /i(z/) = Coo for z/ < z/*. 

2. For /i < Coo and u > u*, u = z/(yu) if and only if /u(z/) = u. 

The reciprocity of /i(z/) and z/(/i) is analogous to the equivalence of the 
grand-canonical and the constant pressure ensembles. Indeed, the pressure 
and the Gibbs energy (per particle) are both obtained as Legendre transforms 
of the free energy, one with respect to the density, the other with respect to 
the volume per particle, 

p{(3, /i) = sup (/ip - /(/5, p)) , g{l3,p) = inf (/(/3, v) + pv) , 
p ^ 

with f{P, v) = vf{(3, v~^) the free energy per particle. Equivalence of ensem- 
bles here means that p{f3, ■) and g{/3, ■) are reciprocal: the Gibbs energy is 
the same as the chemical potential. 

Similarly, p(z^) looks like a Legendre transform oi k Ek with respect 
to A;, while z/(p) looks like a Legendre transform oi Ek/k with respect to l/fc, 
which should be compared with the relations v = 1/p, /(/3,f) = f{f3,p)/p. 



Proof of Lemma A . 1\ 1. The statement for the function p(z/) was proven 



in |5],[T0]. For z/(p), we note that it is the infimum of a family of decreasing, 
affine functions and therefore concave and decreasing. Moreover it is almost 
everywhere differentiable, with derivative —k{fi), the minimizer of Ek — kfi. 
In particular fc(p) > 1, hence z/(/i) is strictly decreasing. 

2. We prove The proof of the converse is similar. Thus let p < Coo 

and z/ = mik{Ek — kfi). Clearly, u < in{k{Ek — kcoo) = v* ■, and for every 
A; G N, 

^ ^ Ek — kii =^ li < — , 

k 

whence p < p(z/). On the other hand, if p < Cqo, then E^ — kp > u* + k{eoo — 
p) — )■ 00 as A; — )■ 00, so there must be a finite k such that u = E^. — kp. It 
follows that p = {Ek — i')/k > piy), whence p = p(z/). If p = Cqo, then 
V = V* and the claim follows from the general inequality p(z/) < Coo- □ 

Lemma A. 2 (Comparison of thresholds). Let 

Pi := mf -, Ui := -pi. 

A:>2 A; - 1 

Then 
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either /xi = Coo and v* = — Coo = ^^i, 
or fii < Coo and z/* < — Cqo < ^^i • 



Proof. Lemma A.l implies the general bounds jji < Coo and z/i > z/*. More- 



over, by definition, v* < Ei — eoo = — Coo and 

ui = sup > lim = -Coo 

k i — K fc-s>oo i — fc 

SO that u* < —Coo < '^i- If in addition /ii = Coo, then i/i = — Coo and for 
all k G N, Ek > {k — l)eoo from which we get u* = mik{Ek — kcoo) > — Coo- 
Since in any case u* < e^o, we get u* = Coo- 

If /ii < Coo, then z/i > — Cqo and there is a p G N such that /xi = Ep/ {p — 
1) < Coo- It follows that 

T^* < Ep- pCoo = ip- - Coo) - Coo < -Coo- □ 

Lemma A. 3 ("Phase" diagram). 1. For /i < /ii, Ej. — k^ has the unique 
minimizer /c(/i) = 1. Similarly, for z/ > z/i, (E'fc — v)/k has the unique 
minimizer fc(z/) = 1. 

2. For /ii < < Coo, every minimizer is finite and larger or equal to 2; 
similarly for v* < v < Vi. 

3. For V < u*, (Ek — v)/k has no finite minimizer. 

Proof. 1. By definition, /i < /ii if and only if for all k > 2, {Ei — = 
< Ek — kn- Thus for ^ < fii, Ek — kfi has the unique minimizer k{^) = 1. 
The statement on {E^ — v)/k is proven in an analogous way. 

2. For /i < Coo, Ek — kfi > k{eoo — /x) — )■ oo as A; — )■ oo, thus {Ek — kfi)km 
reaches its minimum at finite values of k. If k = 1 was a minimizer, we would 
have — /i < Ek — kn for all k > 2, whence /i < yUi. Therefore when /i > /ii, 
every minimizer k{n) is larger or equal to two. The proof for the statement 
on [Ek — T^)/k is similar. 

3. By definition, if u < u*, then v < Ek — kcoo for all A; G N, thus 
[Ek — v)/k> Coo for all k. It follows that in.f k{Ek — i')/k = e^o and there is 
no finite minimizer. □ 



B Free energy at low temperature and low 
density 



Here we give a sketch of the proof of (17). The primary aim is to show that 



Po can be chosen indeed of the order of the preferred ground state density. 
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Po ~ l/a"^. For the sake of completeness, we also make a remark on how 



Eq. (17) should be modified for potentials without attractive tail. 



Potentials with attractive tail Let Z^^{P) be the cluster partition func- 



tion from (29) above. Then [TUl Lemma 3.1] 

N 



zaw,n)< yi n 



The sum is over integers A^i, . . . , Nn G Nq such that kNk = N. The 
integers describe a partition of the N particles into clusters, i.e., groups of 
particles close in space. Using that for suitable c > and all /3 > and 

keN, 

Zl\(3)<exp{-(3E,)expick), 



[ini Lemma 4.3] we deduce that —I3f{(3^p) is upper bounded by the supre- 
mum of 



cp 



p{^^PkEk + (p - '^kpk)ej^ + ^pfc(l - logpfc). 



fceN km fceN 



over all {pk)km ^ [0, oo)^ such that YIT ^Pk < P (think p^ = iVfc/|A|). Next, 
we observe that the mixing entropy can be bounded as J2k Pk^'^siPk/ p) > 
— 2p, for all p > and all admissible (pk), see [TOl Lemma 4.2]. Therefore we 
obtain 

-f3f{f3,p)<{c + 3)p 

- /3inf|(p - ^/cpfe)eoo + J^Pfc(£'fc + /3"Mogp) | ^ A;pfc < p|, 
fceN feeN k 

whence 

/(/3,p) > -(c + 3)rV + pinf (42) 

for all /3 > and all p > 0. 

It remains to obtain an upper bound for the free energy, or a lower bound 
for the partition function. Consider first the case u* > and p > exp(— /3z/*). 



In this case the infimum in Eq. ( 42 ) equals Coo • We lower bound the partition 
function by integrating only over a small neighborhood of the A^-particle 
ground state, and deduce — /(/3,p) > Cqo — C(3~^\ogf3 for suitable C and 
sufficiently large /3. Note that this is possible if p is smaller than the density 
of the ground state, thus p < is sufficient. 
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Next, consider the case p = exp(— /3z/) < exp(— /3z/*). In this case {Ek — 
v)/k has a finite minimizer k = k{v) G N. We lower bound the partition 
function for a cube A = [0, LY and N e /cM particles as follows: we split the 
cube into M small cubes ("cells") with side length of the order ak^^'^ and 
mutual distance R. Here R is of the order of the potential range, and ak^^'^ 
is large enough so that a fc-particle ground state fits into the small cube, as 
in Assumption [2] We can place approximately M = \K\/[ak^/''' + RY small 
cubes in that way. We consider configurations in which particles form clusters 
of size k close to their ground state, such that each cluster fits completely 
into a small cube, and there is at most one cluster per cell. 

We refer the reader to [10] for the details and content ourselves with the 
following remark: the procedure works provided the number M of available 
cells is larger than N/k. This gives the condition 

p < {a + Rk-^/'^)-'^. 

Hence if we choose /3 large enough so that exp(— /3z/*) < (a + R)~'^, the 
condition is certainly fulfilled for every p < exp(— /3z/*). 

Remember that for p > exp(— /3z/*), we are in the first case considered 
above and we only need p < ^ ground state density. Therefore, in the 
end, all we need is the condition p < po with po of the order of 1/a'^. 



Potentials without attractive tail If v has no attractive tail, we might 
have u* = 0, and ground states are not necessarily connected. Set Ef = 
El = and 

El^ := M{U{xi,. . . ,Xk) I {xi,...,Xk) E {R'^Y -R-connected} > Ek- 

The lower bound (42) is still true, but can be improved by replacing Ek by 
El^. In fact, indices k with Ef > Ek can be dropped altogether: if E^^ > Ek, 
then for suitable r > 2, ki + ■ ■ ■ + kr = k, 

Et >Ek = El\ + --- + Ell (43) 

Suppose that for some u > 0, [E^^ — v)/ki > {E'j^ — v)/k for all i. Then 

E > + ^'^^) = - + Et > Et, 



contradicting (43). Thus 



inf = inf{^^ I keN, Ek = Et} 



is the appropriate auxiliary variational problem to be substituted into Eq. (17). 
The density po can be chosen of the order of (a + R)~^^'^. 
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Non-negative potentials When f > 0, the situation becomes particularly 
simple: we have u* = and for all > 0, 



inf 

fceN 



k 



mf — ^- 

km k 



and k{v) = 1 is the unique minimizer. Eq. (17) is replaced by the following: 
For sufficiently low temperature and density p smaller than or of the order 
of l/R^, 

|/(/3,p)-rVlogp| <Cpp-HogP 
and we recognize the free energy [3~^p{\ogp — 1) of an ideal gas. 



C Cluster expansion in the canonical ensem- 
ble 

The virial expansion (|5| can be derived directly with the help of a cluster 
expansion in the canonical ensemble; this was recently done in fTFj. The aim 
of this appendix is to complement Sect, [s] and explain how Eq. (33) (without 
the exact formula for a(m)) is obtained with the approach from |16j . 

The starting point is an expression of the canonical partition function 
as a sum over set partitions {Xi, . . . ,Xr}, r G N, of the particle label set 
{l,...,iV}: 

{Xi,...,Js:,} 

Monomers (|X| = 1) have activity 1, sets with higher cardinality have activity 

\X\ = k>2: Ca(X) = 7^: I W " " " ^^k. 

' ' 7eec(fe) •'^'^ (ii)e7 

with Qc{k) the set of connected graphs with vertices 1, . . . , fc. Note that 
as |A| — )■ oo, for every fixed k and /3, the activity is related to the Mayer 
coefficients as follows: 

The formalism of cluster expansions for polymer partition functions gives 

lAI^ 1 

log Za(/3, AT) = log Li_ + ^ _ Y n{X) Ca(Xi) ■ ■ ■ Ca(X,). 

r>l ' Js:=(X]^,...,Xr) 
conn, XidTj^ 
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Here Fjv is the collection of subsets of {1, . . . , A^} of cardinality at least 2, and 
connectedness and n{X) are defined as follows: With X = {Xi, . . . , Xr) G 
we associate the graph G{X) with vertices 1, . . . ,r and edges i 7^ J, 

Xi n Xj 7^ 0. The polymer X is called connected if the graph of overlaps 
G{X) is, and n{X) e No is the index of i.e., n{X) = n+(X) - n_(X) 

with n±{X) the number of connected subgraphs of G{X) with an even (odd) 
number of edges. 

The A^-dependence in the summation index is slightly inconvenient. We 
remove it by exploiting the invariance with respect to particle relabeling, 

logZ4/).iV) = log^ + f (^)f i ^ nm 

n=2 ^ ^ r=l x=(Xi,...,Xr) 
conn, Xj6r„ 

X Ca(^i) • • • CA{Xr) 1 (U^X, = {1, . . . , n}) . 

In the thermodynamic limit N, |A| 00, for each cluster X = {Xi, . . . , X,.) 
in the sum, 

1 fN\ 1 iV" 

P^(^^JCa(^i)---Ca(^.) ~ ;^|A|i+E^(fc,_i) 

This goes to zero unless 1 + X]i(A;i — 1) — n (note that ">" for every cluster 
X). When this condition is satisfied, the components are necessarily distinct, 
Xi 7^ Xj, and the overlap graph G{X) is necessarily a Husimi graph, i.e., a 
graph whose doubly connected components are complete graphs. Using the 
fact that the index of the complete graph on v vertices is {—l)'"^^{v — 1)!, 
one finds that 

n{x) = {-iy-'f[{v,-iy. 

i=l 

with j the number of doubly connected components of G{X) and Vi, . . . ,Vj 
their respective sizes; thus Vi + • ■ ■ + Vj = r. 

Assuming we can exchange summation and thermodynamic limits, we 
obtain 

00 ^ 

-ISfW, p) = -p{\oip- l) + Y. ^B„(/3) 

» „ r (44) 

n=2 r>2X={Xi,...,Xr} i=l 
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The sum ^'•"^ is over collections of subsets {Xi, . . . , X^} with (Xi, . . . , 
connected, |Xj| > 2 for all z, and such that U^^Xj = {l,...,n} and n = 
1 + E1'(I^.|-1)- Setting 



m 



'3 



|{^ I 1 <^ <r, |X,| =j}| 



we have 



n-l = 5^(|X,|-l) = ^(j-l) 

i=2 



i=l 



i=2 



(45) 



and we obtain the expansion (33) with the information a{m) > 0. Additional 
combinatorial steps would be needed to obtain the formula for the a(rn.)'s, 
but the information that they are strictly positive is all that is needed for 



proofs of Propositions 3.10 and 3.11 
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